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Abstract 



The weak-localization contribution 5cr{B) to the conductivity of a tunnel- 
coupled double-layer electron system is evaluated and its behavior in weak 
magnetic fields B perpendicular or parallel to the layers is examined. In 
a perpendicular field B, 5a{B) increases and stays dependent on tunneling 
as long as the magnetic field is smaller than H/eDrt, where D is the in- 
plane diffusion coefficient and tj the interlayer tunneling time. If Tt is smaller 
than the inelastic scattering time, a parallel magnetic field also leads to a 
considerable increase of the concuctivity starting with a B 2 law and saturating 
at fields higher than h/eZ^/Drt, where Z is the interlayer distance. In the 
limit of coherent tunneling, when t% is comparable to elastic scattering time, 
5o~(B) differs from that of a single-layer system due to ensuing modifications 
of the diffusion coefficient. A possibility to probe the weak-localization effect 
in double-layer systems by the dependence of the conductivity on the gate- 
controlled level splitting is discussed. 

PACS: 73.20.Fz, 73.50.Dn, 73.40.Gk 
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I. INTRODUCTION 



The interference of electronic waves leads to negative corrections to the conductivity of 
electron systems. This effect is known as weak localization. It can be observed at very low 
temperatures, when the inelastic scattering rate is so small that the phase coherence is kept 
over many acts of elastic scattering. Although these quantum corrections are usually small, 
they can be distinguished due to their specific dependence on temperature, magnetic field, 
and frequency of the applied electric field. 

The fundamentals of the theory of weak localization have been developed in Refs. 1-4 and 
a review is given in Ref. 5. Recently considerable attention 6-16 has been drawn to weak- 
localization phenomena in layered systems where characteristic features are caused by a 
dimensionality crossover. These studies have been applied mostly to multilayer systems (su- 
perlattices) although the cases of barrier-separated thin metallic films 9,10 and single-barrier 
structures 12 have also been investigated. In this paper we calculate the weak-localization 
corrections to the in-plane conductivity of two tunnel-coupled two-dimensional (2D) layers, 
the system which is typically formed in double quantum well structures 17 . Such double- 
layer electron systems represent an intermediate case between a single 2D layer and a 
three-dimensional superlattice and their various interesting properties are caused mainly 
by this fact. As concerns the weak localization problem, the tunnel coupling introduces an 
additional cutoff parameter in the diffusion pole, which implies that the weak-localization 
contribution to the conductivity depends on the strength of this coupling. This happens 
when the probability of tunneling is not small in comparison to that of inelastic scattering. 
Another important consequence of the tunnel coupling is that a weak magnetic field B ap- 
plied parallel to the layers leads to a derealization as does a field B applied perpendicular 
to the layers 4 . The physical origin of the effect of parallel field is explained in a similar way: 
this field introduces an additional phase for the electron which moves in one layer, then 
tunnels into another layer, moves there, and finally tunnels back and returns to the initial 
point. The described phenomena are the subject of the present study. 
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The paper is organized as follows. In Sec. II we present the basic formalism for the cal- 
culation of the weak-localization contribution to the conductivity of a double-layer system. 
In Sec. Ill we calculate the magnetoconductivity in a perpendicular magnetic field B and 
in Sec. IV we repeat the calculation for a parallel field B. In Sec. V we study the magne- 
toconductivity in conditions of coherent tunnel coupling, when the tunneling probability is 
comparable to or greater than the probability of elastic scattering. In Sec. VI we discuss the 
results and consider the possibility to probe the weak-localization effects in double quan- 
tum wells by examining the dependence of the conductivity on the gate-controlled energy 
splitting of the lowest two levels of the double quantum well. 



II. FORMALISM 

The Hamiltonian of the double quantum well system is given, in the basis of the left (I) 
and right (r) layer orbitals |Z) and |r), by 

if (x) = Pi [£,(x) + VHx)] + P r [£ r (x) + K(x)] + h (1) 

where x = (x, y) is the in-plane position vector; Ej(x) = (j\ (—ihd/dx. + eA(x, z)) 2 \j) /2m 
(j = /, r) is the matrix element of the kinetic energy operator, A(x, z) the vector-potential, e 
the elementary charge, m the effective mass, and Vj(x) = (j \ V(x,z) \ j) the matrix element 
of the disorder potential. Further, 

h = ya 2 + Ta x , (2) 

is the potential energy matrix of the double quantum well system expressed through the 
energy level splitting A and tunneling matrix element T. Finally, &i are the Pauli matrices 
and Pi — (1 + a z )/2 and P r — (1 — a z )/2 the projection matrices. 

According to the Kubo formula, the zero-temperature dc conductivity is given by 

e 2 h 



27rm 2 L 2 W=R A 



E (-l)Tr/ /^(^ (^^(xOGj ^x)) (3) 
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where G^| A (x, x') are the Green's functions of the electron system with Fermi energy Ep, 
L 2 is the normalization area, 7r a (x) is the kinematic momentum operator, whose matrix 
elements are given as [7r(x)]^ = (j \ — ihd/d~x.+ eA(pc, z) and a and a' are the coordinate 
indices (x, y). The angular brackets (...) denote statistical averaging, Tr denotes the trace, 
and I = 1 for b — b' and I — for b ^ b' . 

The weak- localization contribution 5a aa i to the in-plane conductivity is given by the sum 
of an infinite set of diagrams with two or more maximally crossed impurity lines. Below we 
consider a system of randomly distributed elastic scatterers (impurities) described by the 
correlation function (Vj(x)V^(x')) = W^y(x, x') = Wj8jj>d(x — x'). We obtain 

fcw = 2^TO E E / /^x'4(x)G; ji (x,x 1 )G 3 V(x 2) x') 

bb'=R,A jj'hh 



X7^ i ,(x')G5; /l (x' ) x 1 )G^(x a) x)^ !1 (x 1 ,x j ,), (4) 



where G^,(x,x') = ^[G^(x,x')] jji ) are the averaged one-particle Green's functions, and 
C^- 2 (xi,X2) is the Cooperon, the solution of the Bethe-Salpeter equation, given by 



^iii2( Xl ' X2 ) — ^ji^iuaC*!' x 2)^jii 2 ( Xl ' x 2, 
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+w n £ / dxG^-Cx!, x)G^.( Xl , x)C^ a (x, x 2 ). (5) 

i 

In the following we use Eqs. (f|) and (|5|) to calculate 5a aa > for both directions of the magnetic 
field, perpendicular and parallel to the layers. 

III. PERPENDICULAR MAGNETIC FIELD 

Consider the case when the field B is directed perpendicular to the layers, B = (0, 0, B). 
In the Landau gauge, A = (0, Bx, 0), we can write the Green's function as 

1 



G; ii2 (x 1 ,x 2 ) = — e " - ' <- - r 



X 
n 



i(xi+x 2 )(yi-V2)/2i 2 e -(xi-x 2 ) 2 /4(? 2 

T,L° n ((xx-X2) 2 /2t 2 )G» jlj2 (n), (6) 



where G b j 1 j 2 { n ) is the averaged Green's function in the Landau- level representation, L° are 
the Laguerre polynomials, and I - 



'Ti/eB the magnetic length. From Eqs. (5) and (6) it 
follows that the Cooperon may be written as 



<%(xi,xa) 



<xl+X2){yi . y2)/e 2 C 6^ ( | Xi _ X2 |) ; 



(7) 



where is its translationally invariant part, which is also called Cooperon in the following. 
Introducing C^,{jp) as the Fourier transform of Cj*'(x), we obtain 

C hiM = *>A(P)W» + «fc E / dxeW-*>* 

p'J 

(U - hy/ff + (p' y + Ux/tfY' 2 ) Cf k {p% (8) 

where A^- 2 (p) is the Fourier transform of the translationally invariant part of the product 
(jk - 2 (xi, x 2 )Gj' U2 (xi, x 2 ). In the Landau-level representation it can be expressed as 

Ag(p) = ^ E(" 1 ) n+n ' ^P (-«) Lr 1 ' («) ^ («) ^(rOG^n'), (9) 

tin' 

where u = p 2 £ 2 /2h 2 . Below we consider very weak fields B, I 2 3> (vFjTj) 2 , where VFj and 
Tj are, respectively, the Fermi velocities and elastic scattering times in the layers. In this 



diffusion limit we can expand A*' 



w 
hi 



of Eq. (RJ) in series of hx/£ 2 and hy/£ 2 . From Eq. 



(H) we obtain a set of differential equations 



3 L 



ft 4 r/ 2 A b6 ' (wl r) 2 



2£ 4 dp 2 <9p 



(10) 



In the diffusion limit we can also neglect the B dependence of A^ 2 (jp) , using the expression 

(p) = EgL(p/ 2 - q)^(p/2 + q), (ii) 

q 

where - 2 (p) are the matrix elements of the Green's function 



2m 2 \ 77 r t „i J 2 \r r r ipr/ j 



describing the double-layer system in the absence of the magnetic field. The Fermi energy 
is measured from the center between the two levels of the double quantum well, the elas- 
tic scattering times are introduced as Tj = H 3 /mWj, and t^j are the phenomenologically 
introduced inelastic scattering times. 



In the weak coupling limit 



2{T/h) 2 T i 

< -, (13) 



T t l + (Ar/h) 2 r' 

where r = 2r/r r /(r/ + r r ) is the average elastic scattering time and r t the tunneling time. A 
calculation of A^ j2 (p) in the limit (13) gives 



w>^{ l -%-Mtf- T i) (14) 



and 



A^(P) = A^(p) - J' + JV • (15) 

+ w r )T t 

In Eq. (14) we introduced the 2D diffusion coefficients in the layers Dj = VpjTj/2. A A ?(p) is 
given by Eqs. (14) and (15) as well while Af^(p) and A AA (p) are small and can be neglected. 
A substitution of Eqs. (14) and (15) into Eq. (10) leads to the following equations 

n 2 D lTl d 2 



£ 4 dp 2 



C$ A (P) +(^ + D in ( P -)\ T A C* A {p) - ^C« A ( P ) = w h (16) 



n 2 D r r r d 2 ^ A{p) + (^ + DrTr ( P -)\ T -L) C« A (p) - T ^\p) = 0. (17) 



C^ A (p) and Cfr A (p) are given by Eqs. (16) and (17) with the indices / and r interchanged. 

Equations (16) and (17) can be diagonalized by the substitution Cff A (p) = C[{p) + C l 2 (p) 
and C% A (p) = \iC[(p) + X 2 C l 2 (p). For C k (p) (k =1,2) we obtain 

h 2 d 2 (p^ 2 



£ 4 <9p 2 \h 



Clip) = ^jjA k , (18) 



where 



qf 2 = ( Sl + Sr ) /2 ± S/2, Ai, 2 = 1/2 ± (s, - s r )/2S. (19) 

Here Sj = (Djt^)- 1 + (D^t)' 1 , s 2 = (DiD r T 2 )-\ and S = [(s, - s r ) 2 + 4s t 2 ]V2. Equation 
(18) is analogous to that for the Green's function of a two-dimensional harmonic oscillator 
and its solution is obtained in a straightforward way. The result is 



cr(p) = 2^ r E(- i r^(2^) 



(20) 



_(2/£ 2 )(2n + l) + g 2 (2/£ 2 )(2n + l)+g 2 2 _ 

The expression for C^ r A {p) is given by Eq.(20) with the layer indices interchanged. 

Both terms in Eq. (20) contain diffusion poles modified due to tunneling. In the absence 
of tunneling (uncoupled layers) we obtain the well-known result 4 for each layer. Now we 
calculate the weak-localization contribution 5a aa >. In the diffusion limit Eq. (4) gives 

5a -' = 2^ £ £ E^^,( P )G; 2i ,(p') 

bb'=R,A 33' hh PP' 

x ^(P^>)OlP + P'l)- (21) 

In the limit (|H|) we should retain only the terms with j = f = j\ = j 2 . The weak-localization 
contribution is expressed through C ; f A (p) and C^ T A {p). It gives a positive magnetoconduc- 
tivity Act = 5a(B) — <5er(0) as a sum of two parts 9 

Here 5ao = e 2 /2Tc 2 h, and f(x) = ip(l/2 + 1/x) + hire, where ip(x) is the digamma function. 
In the absence of the magnetic field, the weak-localization contribution, within logarithmic 
accuracy, is given as 



MO) = M In (-^- + 2 ) , (23) 

where we introduced an average inelastic scattering time = 2r^T<p r j '(r v i + r^). It is 
convenient to analyze the results (22) and (23) in the symmetric case described by D\ = 
D r = D, Ti = r r , and t v i = t w . In this case Eq. (23) takes the form 8a(0) = 5ct [In (r/r^) + 
\h.(t/t v + 2r/r t )]. The first term is the same as in the absence of tunneling while the second 
one is tunneling-dependent. In the symmetric case Eq. (22) gives 

M = ; (fib) + ; (flo(l + 2r„/r,)) ' (24) 
where B = h/AeDr^ is the characteristic field. In weak magnetic fields, the magnetocon- 
ductivity increases according to a B 2 law. In stronger fields, when the magnetic length 



i becomes less than y/Dr t , 5a(B) loses its dependence on the tunneling and the system 
behaves like two decoupled layers. In the limit B ^> Bq, B ^> Bolr^/rt) the magne- 
toconductivity is proportional to 25cr lnB. For r t <C t v , an intermediate regime exists, 
B <C B <C Bo(r v /T t ), inwhich the magnetoconductivity is proportional to 5ao\n(B). Fig- 
ure 1 shows the magnetic-field dependence of the magnetoconductivity described by Eq. 
(24) for several values of T^/Tj. 



IV. PARALLEL MAGNETIC FIELD 

Due to the spatial separation of the layers, a field B applied parallel to them can consid- 
erably influence the conductivity of the double-layer system, see, for example, Refs. 18-21. 
This effect occurs as a result of the modification of the electron energy spectrum of the 
coupled layers by the field B. However, under the condition r f ^> r the coherent tunnel 
coupling is suppressed by the elastic scattering and the effect disappears. In contrast, the 
influence of this in-plane field B on the weak-localization part of the conductivity should 
exist as long as r t < t v ; the corresponding characteristic field B\ has to be much smaller. 
An assessment of this influence is given below in the limit (13). 

The parallel field renders the electron Green's function anisotropic. For B = (0, B, 0) we 
choose A = (B(z - z ),0,0) with z = ((/| z \l) + (r\ z |r))/2 and obtain 19 ' 20 



,, , p 2 A„ . ih ( 1 1 \ - ih ( 1 i v 1 



G^(p) 



(25) 



2m 2 ' 2 \ti t v i J 2 \r r V/J 

where A Px = A — uj c Zp x + 5b, ^ c — eB/m is the cyclotron frequency, Z = (r\ z \r) — (l\ z \l) 
is the interlayer distance, and 5b ~ B 2 is a small field-dependent correction to the level 

splitting which can be neglected for small B. 

jjb' 



Using Eq. (p5j) we can calculate A^,(p) from Eq. (pT|). In contrast to the previous 
section, A^'(p) is now anisotropic. In the limit £ 2 ^> ZvpjTj instead of Eq. (14) we have 

A 5 A ( p ) - f 1 - v~ ~ DnKp^PBf+pD/h 2 - T A , (26) 
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where ps = mu: c Z = eBZ, and the upper (lower) sign stands for j = I (j = r). Equation 



(|l5l) remains unchanged. The Bethe-Salpeter equation for the Cooperon now reduces to a 
set of linear algebraic equations 



Op) = wAb>* + E <(p)c^(p) 



(27) 



whose solution is straightforward. A calculation of 5a aa i can be done according to Eq. (pTD, 
with Cj^ 2 (\p + p'|) replaced by C}fj 2 (p + p') found from Eq. (27). Since r t ^> r holds, we 
again put j = j' = j\ = ji an d finally obtain 



Act = Aa sat - 5a V / dx t — j = 



sat 



5cr In 



1 + 



1 + v 
n 



l 



ipr 



(28) 



(29) 



n j \ n j i \ n 

where x = (p/h) 2 and Rj(x) are the fourth-order polynomials defined by 

Ri(x) = x 4 + 2(si + s r — xb)x 3 + [(si + s r ) 2 + 2(s[S r — sf) + 2xb(s t — 2si) + x 2 B ]x 2 

+2[(si + s r )(sis r - s 2 ) + x B (2sis r - sf - s 2 ) + x 2 B si]x + [sis r - s 2 t + x B si} 2 , (30) 

with xb = (2p_e/^) 2 . The result for R r (x) is given by Eq. (30) with the layer indices 
interchanged. Although the field B induces an anisotropy in the Green's functions, the 
localization correction Sa(B) is isotropic because the main contribution in Eq. (21) comes 
from |p| ~ |p'| ~ mvpj. 



Equations (28)-(30) show that in the weak-field region B B\, where B\ = fi/eZ \J Dt^ 
and D = 2DiD r / (Di + D r ) is the average diffusion coefficient, the magnetoconductivity 
follows a B 2 law: 



Aa _ / 2s t eZB\ 2 
5a \ hS J 



si + s r 



2, / si + s r + S 
— m 



(31) 



SlS r — S 2 S Ksi + Sr — S, 

At Tt T v j we obtain Aa/6ao = (B/Bi) 2 . In the opposite limit, a simple expression for the 
magnetoconductivity can be obtained in the symmetric case, Aa/5ao = |(r (/ ,/r t ) 2 (_B/_Bi) 2 ; 
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it shows a substantial suppression of the magnetoconductivity as a result of decreasing 
tunneling probability. 

With the increase of B, the ^-dependence becomes weaker. For B 3> B 1 and B 2 3> 
B 2 Ttp/T t , the integral in Eq. (28) can be neglected and the magnetoconductivity is saturated 
and equal to Acr sai defined by Eq. (29). For r^/r t <C 1 the saturated value goes to zero. 
Recently, a saturation behavior of the weak-localization contribution in parallel magnetic 
fields has been theoretically found for superlattices 15 . 

When the condition r t <C r v holds, one can analytically evaluate the magnetoconductivity 
for B 2 <^ Bfr^/rt, the result is 

Aa = 5a \n[l + (B/B 1 ) 2 }. (32) 

It shows that in the intermediate region B\ <C B 2 <C Bfr^/rt, the magnetoconductivity 
follows a logarithmic law. In Fig. 2 we demonstrate a transition from a i? 2 -behavior at 
weak fields to the saturation at stronger fields for several values of the ratio T t /r v . 

The tunneling rate l/r t can be varied by applying a transverse voltage which changes the 
level splitting A 17 . In the case under consideration Ep Ti/t, a relatively small variation 
of A has no considerable influence on the electron densities in the wells but can dramatically 
modify the tunneling rate. Fig. 3 shows the A- dependence of 5a(B) — <5a(0)| a=o for several 
values of the magnetic field, for both parallel (solid) and perpendicular (dashed lines) fields; 
the thick solid line is the result for B = 0. The curves have a typical resonanse-like behavior 
that is affected differently by the two field orientations. An increase of A always leads to a 
decrease of the conductivity, because in this way the tunneling is suppressed and, therefore, 
the interference is increased. The perpendicular magnetic field tends to smoothen this 
effect because it suppresses the dependence of the weak-localization part of the conductivity 
on tunneling, see Sec. III. On the other hand, the parallel field tends to strengthen the 
tunneling dependence of the weak-localization part: in the saturation regime, cf. Fig. 2, 
this dependence is the strongest. 

Above we calculated the magnetoconductivity in the weak-tunneling regime. Below we 
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evaluate it again in the regime of coherent tunnel coupling. 



V. COHERENT TUNNELING REGIME 

The aim of this section is to study the weak localization effects at sufficiently strong, 
coherent tunnel coupling, described by r t ~ r or even r t <C r. In these conditions one 
always has r t <C t v , which means that there is no competition between the tunneling and 
inelastic scattering processes. As a result, instead of two diffusion poles, see Eq. (20), the 
expression for the Cooperon contains just one pole as for a single-layer system. Hovewer, 
the properties of this pole are somewhat different than those in case of a single 2D layer and 
require a separate consideration. 

Below we assume that Ep is large in comparison with T and A. Introducing a single 
Fermi velocity v F for both layers, we write average diffusion coefficient as D = vpr/2. The 
calculation of the Cooperon according to Eqs. (11), (25), and (27) gives 



(33) 



Wi + w r 

The inelastic scattering time enters the diffusion pole only as an average. On the other 
hand, the Cooperon depends on the asymmetry of the elastic scattering described by the 
parameter /i — (r r — T/)/(r r + 77). In Eq. (33) we have also introduced a A-dependent 
diffusion coefficient according to 

D(A) = — , r(A)= A2 + {k/T)2 ■ (34) 
1 ; l-^A)' 1 ; A 2 + AT 2 + (h/r) 2 { } 

Note that D(A) differs from D only if the elastic scattering is asymmetric. 

If there is no parallel magnetic field (pb = 0), Eq. (33) has the same form as the 
Cooperon of a single-layer system. Therefore, in the diffusion regime Eq. (33) with a 
substitution p 2 — > A%eB(n +1/2) can be directly used to calculate the magnet conductivity 
in a perpendicular magnetic field B. Using Eq. (21), we obtain, in analogy with the single- 
layer case, 



11 



Act = Sa f 



( 



4eD(A)Br, 



) 



(35) 



and, within the logarithmic accuracy, 8a(0) = 5ao In (t/t^,). The only difference between Eq. 
(35) and the well-known expression for a single-layer 2D system is the replacement of D in 
the latter by D(A). The physical meaning of this change is clear and easily understood in 
terms of the resistance resonance phenomenon 22 ' 23 , see also Ref. 17. In the corresponding 
theory 22 ' 23 the factor 1/[1 — fi 2 r(A)j describes the A-dependence of the conductivity and of 
the diffusion coefficient D of the double quantum wells with asymmetric scattering. With 
the increase of |A|, D increases and leads to an increase of the magnetoconductivity Act. 

Now we consider the case of a parallel magnetic field. Using Eqs. (33) and (21), we 
obtain 



The magnetoconductivity described by Eq. (36) does not depend on the scattering asymme- 
try. On the other hand, its dependence on A has the same sign as in case of perpendicular 
field: when A increases, r(A) goes to unity and Act increases. The behavior of the magne- 
toconductivity described by Eqs. (35) and (36) is illustrated in Fig. 4. 

The results for magnetoconductivity described in this section are valid for r t <C t v under 
conditions Dr t eB/h <C 1 for perpendicular field B and Dr t (eZB /h) 2 <C 1 for parallel field 
B. In the limit (13), when r(A) = 1, Eq. (35) is equivalent to Eq. (22) written at r t <C 
and Dr t eB/h <C 1, while Eq. (36) is equivalent to Eq. (32). 



We have investigated the influence of the tunnel coupling betwen two 2D layers on the 
weak-localization-induced magnetoconductivity. This couplng intruduces an extra degree of 
freedom for an electron giving it the possibility to tunnel between the layers and, therefore, 
it reduces the interference effects. As a result, the weak- localization contribution is reduced 



Act = 5a In 1 + r(A)(B/B 1 ) 2 



(36) 



VI. DISCUSSION 



12 



as the tunneling rate overcomes the inelastic scattering rate. In particular, it leads to 
the weakening of the magnetoresistance in the perpendicular magnetic field, see Sec. III. 

The results described in Sec. Ill are similar to those obtained in Ref. 9 for two thin 
metallic films divided by a barrier if one equals Tt to the tunneling times T\i and r%i of Ref. 
9. This is not surprising. When the inelastic scattering lengths are small in comparison 
to the film widths, the diffusion of the electrons in such a system proceeds in the same 
way as in coupled 2D layers. A characteristic feature of the double quantum well case, 
as opposed to the case of thin films, is the strong resonance dependence of the tunneling 
time Tt, cf. Eq. (|T^), on the energy-level splitting A. The latter can be easily controlled 
by external gates enclosing the double quantum wells. This opens a new way to examine 
the weak-localization phenomena by measuring the A-dependence of the conductivity, see, 
for example, Fig. 3, in addition to the B dependence. The conductivity, due to its weak- 
localization part, is maximal when the tunneling resonance condition A =0 holds (this can 
be called "delocalisation resonance"), while an increase of A increases the interference and 
leads to a decrease of the conductivity down to its value for uncoupled layers. On the other 
hand, if the scattering times in the layers are different, the main part of the conductivity 
will show the opposite effect (the resistance resonance 17 ). Although the relative size of 
the resistance resonance would be small due to the parameter r/r t , it can obscure the 
derealization resonance. Therefore, the samples with symmetric scattering are preferrable 
for probing the interference corrections. 

In Sec. IV we studied the magnetoconcuctivity in a parallel magnetic field. Since this 
magnetoconcuctivity can exist only in the presence of the tunnel coupling, it becomes more 
pronounced as the tunneling rate overcomes the inelastic scattering rate. Although at very 
small B both parallel and perpendicular fields give £? 2 -corrections to the magnetoconduc- 
tivity, at higher fields the behavior becomes completely different. When the perpendicular 
field is larger than both Bq and BoT^/rt, the weak-localization contribution logarithmically 
increases with B and does not depend on tunnneling. When the parallel field is larger than 
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both Bi and Bi^Jr^/rt, 5cr(B) becomes saturated, i.e., it loses its dependence on B but not 
that on tunneling. Since (vf/Z)^/t^t 3> 1 holds, the characteristic parallel field B\ is much 
larger than the characteristic perpendicular field B . 

With an increase of A both 5a(B) and Act in a parallel magnetic field decrease. However, 
in the case of coherent tunnel coupling (Sec. V) the behavior of the magnetoconductivity Aa 
versus A is opposite: an increase of A increases Aa. A transition between these two kinds 
of behavior can be explained as follows. When the tunnel coupling becomes so strong that 
DT t (eBZ/h) 2 «1, a change of r t with A no longer modifies the magnetic-field suppression 
of the interference, see Sec. IV. With a further increase of the tunneling rate, when the latter 
overcomes the elastic scattering rate, coherent coupling effects become important. Instead 
of the states localized in the wells, one has a pair of hybridized states, which are well-defined 
in the limit of 2T 3> h/r. At A = the electron density in both these states is distributed 
equally among the wells, and the parallel magnetic field cannot reduce the interference in 
the way described in the Introduction. The magnetoconductivity is zero, which is given also 
by Eq. (36) at 2T > /t/r and A = 0. At A ^ the symmetry is broken. The two states are 
mostly localized in different wells and the parallel magnetic field can reduce the interference 
because the transitions of electrons between these states are possible. In conclusion, the 
magnetoconductivity increases with an increase of A. 

Let us briefly discuss the main approximations used in this paper. Since we neglected the 
electron-electron interaction, which also leads to quantum corrections to the conductivity, 
the theory is applicable at temperatures considerably lower than h/r^ 24 . We also neglected 
the possibility of spin-flip scattering, assuming that the relevant scattering time is larger 
than Tp. The assumption of the 5-correllated scattering does not lead to principal changes. 
For example, if one takes into account the interlayer correlation of the scattering amplitudes, 
one has to replace only 1/r in Eq. (13) by 1/r — 1/r^, where r/ r = h 3 /mwi r . 

Throughout the paper we calculated the total, i.e. summed over both layers, in-plane 
conductivity of the double-layer system. Experimentally it corresponds to the double quan- 
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turn wells with common contacts to the layers. The case of independently contacted layers 
does not require a separate consideration for the following reason. Since the system is 
assumed macroscopic, its in-plane size should substantially exceed the inelastic scattering 
length vft^. On the other hand, the tunnel coupling manifests itself in weak-localization 
phenomena under the condition > r t . In such characteristic current redistribu- 

tion length 25 l s ~ \jvpTTt is small in comparison to the size of the system, which means that 
for any case of independent contacting the double-layer system would behave as a system 
with common contacts. 
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FIGURES 

FIG. 1. Magnetoconductivity of the double-layer system in a perpendicular magnetic field 
B for symmetric conditions. The curves are marked with the values of the ratio n/r^ and 
B = h/AeDr^ 

FIG. 2. Magnetoconductivity of the double-layer system in a parallel magnetic field B for 
symmetric conditions. The curves are marked as in Fig. 1 and B\ = h/eZ^/Dr^ 

FIG. 3. Relative weak-localization correction to the conductivity of the double-layer system 
as a function of the level splitting A for several values of the parallel (solid) and perpendicular 
(dashed) magnetic field. It is assumed that Tt/r^ = 0.1 at A = 0. The numbers next to the curves 
show the values of B/Bq (perpendicular field) and (B/Bi) 2 (parallel field). 

FIG. 4. Magnetoconductivity of the double-layer system in a coherent tunneling regime for 
2Tr/h = 4 and ji = 0.7. The solid and dashed lines correspond to parallel and perpendicular field, 
respectively. The curves are marked with the values of the ratio A/2T. 
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Fig. 1 O.E.Raichev and P.Vasilopoulos, "Weak localization..." 




Fig. 2 O.E.Raichev and P.Vasilopoulos "Weak localization..." 
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Fig.3 O.E.Raichev and P.Vasilopoulos "Weak localization..." 




Fig. 4 O.E.Raichev and P.Vasilopoulos "Weak localization..." 



